To appear in ApJ 



00 
O 

o 

(N 

a 

S 1 

Oh- 

6 

o3 



The most massive core collapse supernova progenitors 

R. Waldman 

Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 

waldman@cc . huj i . ac . il 

ABSTRACT 

The discovery of the extremely luminous supernova SN 2006gy, possibly inter- 
preted as a pair instability supernova, renewed the interest in very massive stars. 
We explore the evolution of these objects, which end their life as pair instability 
supernovae or as core collapse supernovae with relatively massive iron cores, up 
to about 3 Mq. 
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Subject headings: stars: evolution, (stars:) supernovae: general 



1. Introduction 

The interest in the evolution of very massive stars (VMS), with masses > 100 M , 
has recen tly been revived by the discovery of SN 2006gy - the most luminous supernova ever 
recorded (jOfek et al.ll2007l : ISmith et al.ll2007l ). This object, having a luminosity of ~ 10 times 



that of a t ypical core-collapse SN (CCSN), is probably the first evidence of a pair instability 



SN (PISN) I Woosley et al.l (120071 ). PISN are massive stellar objects, whose evolutionary path 
brings their center into a region in thermodynamical phase space (p < 10 6 ,T > 10 9 ), where 
thermal energy is converted into the production of electron-positron pairs, thus resulting in 
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ever the overall interest in this topic has been relatively small, mainly due to lack of obser- 
vational data. 

It was originally believed that stars massive enough to produce PISN could only be 
found among population III stars with close to zero metallicity (Z < 10~ 4 ), and hence 
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only at very high redshift (z > 15). More recently IScannapieco et al.l (120051 ) discussed the 



detectability of PISN at redshift of z < 6, arguing that metal enrichment is a local process, 



theref ore metal-free star-forming pockets may be found at such low redshifts. lLanger et al. 



(120071 ) introduced the effect of rotation into studying this question concluding that PISN 
could be produced by slow rotators of met allicity Z < Z^/3 at a rate of one in every 
1000 SN in the local universe. Furthermore, ISmith et al.l (120071 ) point out, that mass loss 
rates in the local universe might be much lower than previously thought, so that massive 
sta rs might be left with e nough mass to become PISN. This conclusion is also supported 
by lYungelson et al.l (120081 ) who extensively discuss the mass loss rates and fates of VMS. 
It is interesting to note , that SN 2006gy t o ok pl ace in the nearby Universe. Following the 
discovery of SN 2006gy, lUmeda &: Nomotol (120081) addressed t he que stion of how much 56 Ni 



can be produced in massive CCSN, while iHeger fc Woosleyl (120081 ) computed the detailed 
nucleosynthesis in these SNe. 

The interest in VMS is further motivated by the discovery of Ultraluminous X-ray 
Sources (ULX), which can be interpreted as mass-accreting intermediate mass black holes 
(IMBH) with mass ~ (10 2 — 10 5 ) M & . One of the possible scenarios for IMBH formation is by 



VMS formed by stellar mergers in comp act globular clusters (see e.g . lYungelson et al.l 12008 



and references therein). In this context, iNakazato et al.l (120061 . 120071 ) studied the collapse of 
massive iron cores with M > 3M . In their first paper they treat the fate of stars of mass 
> 300 M & which reach the photodisintegration temperature (~ 6 x 10 9 i^) after undergoing 
pair instability. The entropy per baryon of these models at photodisintegration is s > 16ks 
compared with the classical core-collapse SN with s ~ Iks- In the second paper they aim to 
bridge this entropy gap, corresponding to core masses of (3 — 30) M & but claim that there is 
a lack of systematic progenitor models for this range, hence they use synthetic initial models 
for their calculations. 

In this work we focus mostly on the mass range M < 80 M (He core mass M# e < 
36 M Q ) immediately below the range which enters the pair instability region, and present a 
systematic picture of the resulting CCSN progenitors. 



2. Method 



Since the mass loss r ates of stars in this range are highly uncertain, (see e.g. discussion 
by lYungelson et al.l 12008!). we avoid dealing with this question by following the example 
of IHeger &: Woosleyl (120021 ). and modeling the evolution of helium cores. Our helium core 
initial models are homogeneous polytropes composed entirely of helium and metals, with 
metallicity Z w 0.015, in the mass range (8 — 160) M & . The models were then evolved to 
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the helium zero age main sequence. In the following we will refer to these models as "HeiV" 
where N is the mass of the model. For comparison we evolved also a few models of regular 
hydrogen stars, beginning from the zero-age main sequence (ZAMS). We will refer to these 
models as "MiV" where N is the mass of the model. All our models have no mass-loss. We 
argue that as long as the mass loss rate is not so high that it will cut into the He-core, 
the evolution after the main-sequence phase will be virtually independent of the fate of the 
hydrogen-rich envelope. We followed the evolution of each model until the star is either 
completely disrupted (for the PISN case) or Fe begins to photo-disintegrate (for the CCSN 
case) . 

We followed the evolution using the Lagrangian one dimensional Tycho evolutionary 
code version 6.92 (with some mod ifications), publicly available on the web (the code is 
described in lYoung fe Arnettl 120051 ). Convection is treated using the well known mixing 
length theory (MLT) with the Ledoux criterion. In the MLT f ormulation of Tycho, th e value 
of the mixing length parameter fit to the Sun is oimlt ~ 2.1 (lYoung fe Arnettl 120051 ) . so we 
used a value of oimlt = 2 in our calculations. The nuclear reaction rates used by TYCHO are 
taken from the NON-SMOKER database as described in lRauscher fe Thielemann! (120001 ) . 



The evolution is generally followed using the code's hydrostatic mode. The pulsational 
pair instability models are treated as follows. When hydrodynamic instability is encountered, 
the code is switched to the hydrodynamic mode, and mass ejection is accounted for by 
removing outer zones having supersonic velocity in excess of the escape velocity. After mass 
ejection has died out, and the stellar core is already in contraction, the code is switched back 
to the hydrostatic mode to follow the interpulse period. 



3. Results 

The He-core models we computed can be divided into four categories, according to their 
final fate, as can be seen in the central density and temperature plot (Fig. [TJ): 

1. CCSN - Models that reach core collapse (i.e. Fe photo-disintegration) conditions with- 
out entering the region of pair instability. This is the fate of He-cores with mass 
M < 36 Mq, as can be seen for the models He8 and He36 in Fig. [TJ 

2. Pulsational PISN (PPI)- Models that reach pair instability, collapse and bounce due 
to the energy released by nuclear reactions, but the energy released is insufficient to 
disrupt the entire star, thus a fraction of the star's mass is emitted, and the star 
collapses back. This may happen several times, until the star has no more material to 
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burn and reverse the collapse, and core collapse conditions are reached. This occurs 
for models with He-core mass in the range 36 < M < 54 M , e.g. model He48 in Fig. 

m 

3. PISN - Models that reach pair instability, collapse, and the energy released by nuclear 
reactions is high enough to disrupt the entire star. This occurs for models with He-core 
mass in the range 54 < M < 130 M , e.g. model He80 in Fig. [TJ 

4. Pair instability core collapse (PICC) - Models that reach pair instability, but the energy 
released is too low to reverse the collapse, and the star continues collapsing into the 
photodisintegration regime. This occurs for models with He-core mass in the range 
M > 130 M Q , e.g. model Hel60 in Fig. [TJ 
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Fig. 1. — Evolution of the central density and temperature. Each line is labeled "M" for 
stellar models and "He" for He-core models, followed by the mass of the model. The figure 
is divided into two panels for clarity - the left panel shows models that reach CC without 
reaching pair instability, the right panel shows models reaching pair instability, subsequently 
experiencing pulsations (He48), complete disruption (He80), or direct collapse (Hel60). 
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The properties of our models at pre-SN are summarized in table HJ 



Table 1. Properties of the pre-SN models. 
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Note. — The columns represent for each model the total mass (M), He-core mass (M[j e ), CO-core mass (Mco)> Si-core mass 
(Msi), Fe-core mass (Mp e ), 12 C and 20 Ne mass fraction at the end of core He burning (Xci2 and Xjv e 2o)i central density 
(pc), temperature (T c ), entropy per baryon (5 C ), electron mole fraction (Ye), and binding energy (Eti n ). Masses are given in 
Mq, density in g cm -3 , temperature in 10 9 K, and energy in 10 51 erg. For the models reaching core collapse the data are given 



when central temperature reaches 7 X 10 9 , while for the models that disrupt after reaching pair instability the data are given at 
onset of instability. 
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Fig. [2] shows the density structure of the pre-SN, at the moment when the central 
temperature reaches 7 x 10 9 K. The two extreme models He8 and He36 are shown, as well 
as M80 which has a He-core mass similar to the He36 model, and M20 - a typical CCSN 
progenitor. The composition of the same models is shown in Fig. [3J 
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Fig. 2. — Pre-SN density structure. Each line is labeled "M" for stellar models and "He" for 
He-core models, followed by the mass of the model. (In the color version the solid part of each 
line designates the Fe-group core, the dotted part - the Si-group core, and the dash-dotted 
- the rest of the model.) 
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Fig. 3.— Pre-SN composition of models M20, He8, M80 and He36. "Si" and "Fe" stand for 
the total of Si- and Fe-group elements respectively. 
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The size of the Fe-, Si-, and CO-cores of our pre-SN models is shown in Fig. H] plotted 
against the size of the He-core. A scaled-up view of the size of the Fe-core (defined as 
the mass coordinate where the electron mole fraction Y e < 0.49) together with the central 
entropy per baryon for the same models is shown in Fig. Note that the size of the Fe-core 
is slightly non-monotonic. The central entropy, is monotonic with mass, but slightly differs 
between He-core and stellar models. 

From the above results it is notable that the He-core models behave similarly to the 
stellar models (compare e.g. models He36 and M80 which has a He-core mass of ~ 36 M Q ), 
however some differences still exist (e.g. the central entropy in Fig. [5]). 
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Fig. 4. — Mass of the Fe-core (squares), Si-core (triangles), and CO-core (circles) for the 
computed models. Filled shapes designate He-core models, open shapes - stellar models. 
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Fig. 5. — Mass of the Fe-core (squares) and central entropy per baryon (circles) for the 
computed models. Filled shapes designate He-core models, open shapes - stellar models. 
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Conclusions 



Our results are in gene r al agreement with previou sly published results (e.g. lHeger fe Woosley 



20021 ; IWoosley et al.l 120071 ; lUmeda fc Nomotol 120081 ). We focused on the heaviest models 



which do not encounter pair instability (CCSN) in the range (He8 - He36). The outstanding 
novel features of these models are: 

1. Relatively large Fe-cores, up to about 3M Q , and a large amount (up to about 10 M Q ) 
of Si-group elements. 

2. Comparatively low central density and high central entropy. 

3. A comparatively shallow density profile. 

These differences might have a considerable impact on the behavior of these models 
during core collapse and on the outcome of the explosion, a question which we hope to 
address in the future. 

Similar features are encountered for the lower mass part of the pulsational pair instability 
models (He38 - He50). However, due to the numerical complexity of following the pulsations 
and the potential sensitivity of the results to the numerical treatment of convection, mass 
loss etc., further work is needed to ascertain the validity of the results for the pulsational 
models. 



I would like to thank Zalman Barkat for many hours of fruitful discussion. I would like 
to acknowledge David Arnett for providing his TYCHO code for public use. 
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